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Abstract 

A leading-twist factorization formula is derived for the longitudinal structure function 
in the re — > 1 limit of deeply inelastic scattering. This is achieved by defining a new jet 
function which is gauge independent and probes the transverse momentum of the struck 
parton in the target. In moment space, terms of order In'^ N/N, which are the leading ones 
for Fl, are shown to be resummable through the cusp anomalous dimension and the 
anomalous dimension 7j/ of the new jet function. This anomalous dimension is computed 
to 0{as)- The suggested factorization for reproduces the fixed order results known 
to 0{a1). The general ideas for resumming the terms of order In'^ N/N in moment space 
may be extended to the other structure functions and to other inclusive processes near 
the elastic limit. 
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1 Introduction 



Typical applications of QCD perturbation theory involve processes with two mass scales, the 
momentum transfer Q and the scale of strong interactions Aqcd with ^ ^qcd- Within this 
large family of hard processes a particular class is of great theoretical and phenomenological 
importance. These are the quasi-elastic processes, characterized by a third scale M^, the 
invariant mass of the final states, with ^ M^. Here perturbation theory can be used to 
study the transition from short- distance dynamics (M^ S> Aqqj-,) to long-distance dynamics 
as the elastic limit is approached (M^ Aqqj-,). Deeply inelastic scattering (DIS) structure 
functions at large Bjorken-x, Drell-Yan production near the partonic threshold and thrust in e^- 
e~ annihilation near the two-jet limit are all characteristic examples of quasi-elastic processes. 
In this paper we focus on the large-x DIS structure functions for definiteness, but we expect 
our discussion to be generalizable to all quasi-elastic processes. 

Near the elastic phase-space region x — >■ 1 of DIS the coefficient functions receive corrections 
in 1—x. For F2 the leading power in x corrections are terms of the form a"[ln"*(l — x) / (1 — + , 
with m < 2n — 1. The formalism of resumming this class of corrections has been developed 
in Refs. 121 E] and recently reviewed in Ref. j3]. A key step is the factorization of the 
perturbation theory diagrams at large x in terms of a hard amplitude, a universal soft function, 
and a jet function. We call this factorization "light-cone expansion", distinct from the Wilson 
operator product expansion. The soft function and the jet functions contain singular pieces 
like ln'^(l — x)/{l — a;)+ and 1/(1 — x)^ respectively, which are resummed. All these corrections 
correspond to 0{a'^\n^ N) terms in moment space. 

The question we consider here is whether a formalism beyond the light-cone expansion can 
be constructed to resum terms that do not correspond to power singularities in x and are 
still leading power in Q. A class of such terms are the simple logarithmic corrections of order 

ln™'(l — x), which in moment space generate 0{a'^\n'' N/N) power suppressed contributions. 
The cleanest context where such corrections are encountered is in the longitudinal structure 
function Fl- Fl, as distinct from F2, does not contain any singular distributions in x. Its leading 
power in x corrections are all logarithmic and they start at order a^. A phenomenological 
motivation for studying the resummation of Fl near the elastic limit comes from considering 
the ratio R of the longitudinal and the transverse cross sections : 



where Fj, i = 2,L denote the structure functions, M is the target mass and x = —q^/{2p ■ q). 
The longitudinal structure function Fi vanishes due to helicity conservation in the naive parton 
model and thus R provides us with a clean test of perturbative QCD. There already exist some 
experimental data [HI on this quantity and more is expected in the near future. Leading 0{as) 
and next to leading 0{a1) perturbative corrections to F^ have been calculated [Z]-|n] including 
the target mass corrections. Comparison with the data indicates that the theoretical predictions 
lie systematically below the experimental points, especially in the large x region. This may be 
indicative of rather large higher twist corrections ^01- In order to shed more light on this issue 
it is worth studying the effects of resummation at the leading twist level. 

The resummability of quasi-elastic corrections that do not correspond to singular distri- 
butions was considered in Ref. [Sj. There it was argued that final state power suppressed 
corrections of order a"(l — xY \n^{\ — x), with p > 0, can be identified and resummed through 
an expansion of the final state interactions in terms of higher dimensional bilocal operators. 



R 



Fl 
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The expectation values of these operators define classes of new jet functions that generalize the 
jet function of the light-cone expansion. In this new expansion operators of different dimen- 
sionality all contribute to the leading power in Q (leading twist) after integration over the final 
state phase-space. However, higher dimensional operators produce terms suppressed by higher 
powers of 1 — X near the elastic limit x ^ 1. In Ref . jH] rules were given for identifying indepen- 
dent sets of these higher dimensional operators. Once they are constructed, their anomalous 
dimensions can be used to resum the power suppressed logarithmic corrections, just like in the 
case of the light cone expansion. 

Within the context of the above formal construction, in this work we consider the longi- 
tudinal structure function in the elastic limit as a specific example of a quasi-elastic process. 
We shall show that by analyzing the perturbative diagrams in the appropriate infrared lim- 
its we, indeed, obtain the operator expansion suggested in Ref. jH]. We accomplish this not 
by dimensional arguments but through a diagrammatic analysis by explicitly considering the 
subleading terms in the eikonal approximation for the final state. The result is a factorization 
formula for Fl that is used to resum in moment space all the 0{hi^ N/N) corrections. We find 
it remarkable that such power suppressed corrections can be treated systematically using the 
standard methods for analyzing the infared region of perturbation theory. Simple logarithms 
of 1 — X come from the pinch singular surfaces of Feynman diagrams. This observation, we 
believe has not been considered in the literature so far. 

The general ideas for resumming terms of order In'^ N/N in moment space can be extended 
to other structure functions and to other inclusive processes near the elastic limit. For example, 
it has been noticed ^1] in the context of Higgs production that such corrections can be phe- 
nomenologically important. They may also be of relevance in renormalon phenomenology fT^ . 
Thus, the analysis of the resummation of the O {In'' N/N) corrections is easiest done for the 
case of the longitudinal structure function, but as discussed, the applications go well beyond 
that to phenomenologically more interesting situations which we hope to address in the future. 

The structure of this paper is as follows. In section 2 we set notations and review known 
fixed order results. In section 3 we briefiy summarize the light-cone expansion for F2, which is 
used as a basis for what follows. Then we go on to show how the new jet function, characteristic 
of Fl, arises from the infrared singularities of Feynman diagrams, thus explicitly constructing 
the relevant jet operator of Ref. Here we also study the renormalization group properties of 
the functions entering the Fl factorization formula. In section 4 we present an explicit one-loop 
calculation of the anomalous dimension of the new jet function, which is subsequently used in 
the Fl resummation formula. The resummed coefficient functions are found in agreement with 
the fixed order results. We summarize our conclusions in the final section. Certain useful 
formulae are listed in the Appendix. 

2 Notations and known results 

We consider the DIS structure functions F2,Fl of a massless quark of electric charge Qf = 1 
and with collinear divergences dimensionally regularized in D = 4 — 2e dimensions. Their 
perturbative expansion is 

F,(x,Q^e) = ^ p FW(x,g^e), r = 2, L . (1) 
n=i ^ vr / 
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The DIS quark tensor for an electromagnetic probe is defined as 



W,,{va) = ^J d^'ye-'^y {p\jliO) Uy)\p) = ^J dPSn{M},ML) , (2) 

where PSn is the Lorenz invariant phase space for n particles in the final state and M.l,r are 
the left and right production amplitudes. Structure functions are obtained via the projections 
F, = PI"" W^^, for r = 2, L, given at D = 4 by 

Pr = ^P'P-", Pr = -V'" + ^Pr- (3) 

The normalization in Eq. Q is chosen so that ^2°'* = 6{1 — x). 

The light cone is defined by the two dimensionless vectors f and w'^ normalized as f ■ = 1 
and we fix the v'^ (+) direction to be parallel to the incoming quark momentum p. Then the 
(— ) direction is set by the vector q = xp + q. We denote by pf the total final state momentum, 
Pf = p + q. The light cone (Sudakov) decomposition of the external momenta is 

p/^ = -^^^ p =(l_a;)-^^;'^ + -^^;^ q=-^v'^, (4) 
y/2 ^ ^ ' ^ ^/2x ^X ^ ' 

and for every momentum vector 

k^ = {v- k)v^ + {v ■ k)v^ + k^^ , d^k = d{v ■ k) d{v ■ k) d^-^k^. (5) 

The two massless particle phase space in this basis and for azimuthically symmetric integrands 
becomes 

with the Sudakov variable /32 defined as P2 = ■ k2)/{v ■ q). The three massless particle phase 
space in terms of Sudakov variables is given in the appendix. 

At high momentum transfer, ^ ^qcd) ^"^^ fixed x all structure functions satisfy the 
Wilson operator product expansion. At leading twist this can be conveniently written as a 
factorization theorem in terms of the convolution symbol 



{h®f){x)= / dxi / dx2S{x — X1X2) h{xi) f{x2)- (7) 
Jo Jo 

Specifically for Fl we have 

Fl{x, Q\ e) = Cl{x, Q''lii\ a,{i?)) ® f{x, e, ^.(/i^)) + 0(A^cd/Q'), (8) 

with fl the factorization scale and / the quark distribution function of the quark target. The 
coefficient functions are determined given a factorization scheme, which will be the MS scheme 
for this paper. The coefficient function Cl has been calculated to 0{a'j.) and the results are 

lainiiini 

C^\x,Qyfi') = CFX, (9) 
■1 ■ T . n2 



X] 



ln% + Cf (10) 
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with the well known functions 



qq 



and 



(x,l) 



2 

+ 



-Ci ln"(l-x) 



Pi 



2C(2) 



(C(2) - 1)CfCa - -PiCf 



Cp\n{l-x) 
ln(l — x) + reg. 



:i2) 



In the above expression only terms that are singular in the x — > 1 limit have been retained. It 
is these logarithms of 1 — a; that are the focus of the present study. 



3 The factorization formula for Fl 

We shall show that Fi satisfies a novel factorization formula in the x — > 1 limit which can be 
used to sum the singular ln'^(l — x) contributions in the coefficient function, (see Eq. (fT2|l ). 
This factorization formula involves a new jet function that contains these leading corrections 
to Fl in the large x limit, which in moment space are of 0{\n'' N/N), N being the moment 
variable. Such contributions are also present in F2, although they are subleading to the ones 
that are of 0{\n'' N) in moment space. 



3.1 The new jet function 

Before presenting the factorization formula for F^ in the x —>■ 1 limit it is worth reviewing 
the light-cone expansion for F2, since the latter will serve as a prototype for the former. The 
leading power in x factorization of F2 can be cast in the following form 

F2(x,Q2) = \H2{Q^)\^ dx' dw J[{x' -x-w)Q^) V{w)4>{x') 

= \H2{Q^)\^J^V^4>. (1) 

H2[Q'^) is the hard scattering function on either side of the final state cut. |if2(Q^)P is by 
construction UV dominated, i.e. all its internal lines are far from their mass shell. V is the 
universal soft radiation function defined as 

V{w) = |d|/-e-™^"P-^(0|$t(0,-oo)<l>,(0,y-t;)<l>,(r^,-cx))|0) 

= FtL'I, (0|$|;(0, -00) $,(0,^-1;) $,(y-i;,-oo)|0). (2) 

FT^^p stands for the Fourier transformation along the light-cone direction —v. The Wilson 
line operator starting in the far light-cone past and defined along the direction is given by 

^v{x + tv, —00) = P exp (^—igs J ds v'^ Af^{x + sv)^ . (3) 

Its perturbative expansion generates the Feynman rules for the eikonal approximation, (j) is the 
light-cone parton distribution function. Its operator definition, for the scattering of quarks, is 

0(x')=FTg(p|^(O) ^$,(0,i/-)^(2/)b)®\/-S (4) 
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This distribution function is not identical to / used in Eq. (jH}, but they both share the feature 
that the only singularities they contain come from coUinear propagation along the initial [v) 
light cone direction. The role of in Eq. (jH) is to remove soft contributions, which would 
otherwise be double-counted in both and V. 

Finally, J is the current jet function that describes the stream of almost coUinear particles 
with invariant mass {x' — x — w)Q'^ moving along the v light-cone direction. J is defined as a 
gauge invariant and soft subtracted fermion propagator 

J ((1 - z)Q') = FT^Xp (OI-^-UO, -oo) ^(0) i^iy) -oo)|0) V'', (5) 

where FT^+^p denotes the four-dimensional Fourier transformation with respect to the total jet 
momentum q + zp. 

We note that the factorization theorem of Eq. is written in terms of UV renormalized 
effective non-local operators. When one analyzes the contributions to the characteristic func- 
tions 1-^2^, V, J and coming from classes of graphs of the full theory, it is often the case 
that the assignment to soft and jet functions of these contributions is gauge dependent with 
gauge independence occuring only for their convolution, i.e. the final result. Here, however, the 
characteristic functions are defined in terms of gauge invariant operators, thus we can classify 
singularities and enhancements as soft or collinear in an unambiguous way. The cost for this 
is the need to perform infrared subtractions, as denoted for instance by the denominator V in 
Eq. (0) for the jet function. 

In Ref. jS] it was argued that Fl satisfies the following factorization formula: 

Fi(a;,Q2,e) = |ifi(Q2)|2j'®\/®0. (6) 

The factors V and 4> are the same as in the F2 case, see Eqs. (0), ©• The new jet function J', 
is defined as the expectation value of a bi-local operator that probes the transverse momentum 
of the outgoing struck quark. Explicitly, 

f((l-z)Q') = (^-^^^ FT^XpmliO^-^)PMO)P±m^viy,-oo)\0)®V-\ (7) 

The conventional normalization factors of Eqs. (jH)) have been included in this definition. 
The operator defining J' is the only gauge invariant bi-local operator that is non-vanishing by 
the massless quark equations of motion or the projection of Eq. Q, and although it is of 
higher dimension than the corresponding one that defines J, Eq. (0), it contributes to leading 
twist. All the terms that it generates are of order (1 — x)^, i.e. only logarithms. It is therefore 
the appropriate operator for resumming the leading logarithmic quasi-elastic corrections in Fl. 
In Ref. the new jet operator was identified by requiring it satisfied all the above properties. 
In the next subsection we show that it actually comes out naturally from the infrared properties 
of the Feynman diagrams that contribute to F^ at leading twist. 

The jet function J' is depicted in fig. 1(a), and the factorization expression for F^ is depicted 
in fig. 1(b). We emphasize that the factor V (p is common to every DIS observable, and that 
the only difference between F2 and Fl at their respective leading power in x level is that they are 
different Lorentz projections of the same tensor W^uip, q), hence they differ in the jet functions 
entering their respective factorization formulae. 
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3.2 Operators beyond the eikonal approximation 

In order to prove the factorization formula in Eq. © in this subsection we follow the constructive 
approach. We analyze the perturbation theory diagrams and show how the operator in Eq. ((7j) 
emerges once one goes beyond the eikonal approximation for Fl. Infrared singularities and 
ensuing corrections arise from characterisitc regions in momentum space called pinch surfaces. 
These are the same for all observables, therefore the classification of internal lines in an Fl 
diagram are the same as in any other observable, i.e. we have collinear, soft and ultraviolet 
lines. The factorization of the soft lines from the incoming f-jet is as in F2. We, therefore, 
concentrate on the factorization of v-coUinear lines from the initial state w-jet since this is where 
the difference between Fl and F2 lies. 

Consider the configuration shown in fig. 2, in which the depicted gluon belongs to the w-jet. 
Clearly the fermion line carrying momentum (p — I) is far off-shell and belongs to the hard 
scattering subdiagram in the notation of reduced diagrams. If, however, we were to factorize 
the collinear gluon as usual by eikonalizing the fermion propagator then this, in general leading, 
IR contribution would be projected out in F^ because of the contraction with at the current 
vertex. Note that also in the case where the gluon is soft, the leading eikonal contribution is 
similarly projected out. Therefore, we must go beyond the eikonal approximation. 

The graph in fig. 2 is proportional to 

M- - 'J'^, f = * (^^ - ^ /) + A (8) 

We define the transverse projector lo^ = 5"^ and decompose every 7 matrix in light cone 
components as 

= v^ /j + v^ p+ (tu^)^. (9) 

We also define the effective vertex 

O"(0=^7.((^7)'-(/^7)^)- (10) 
It enjoys the following properties 

O\l)h = 0, O\l)v, = 0, O\l)vx = -gsilcoi)^^, 0\l)Oxin = -2gl (11) 

and 

0\l)...0,{n=g',{u;^)\..{u;^),. (12) 
Then the factor in Eq. (jHJ becomes 

0\l) - l-0\l) p /J + g^v' /, + g,J—{luj^){u;^)' (13) 
2 2v ■ I 

If only one gluon from the "W-jet is connected to the hard scattering, as in fig. 2, then only the 
first term in (fT!?|l survives the massless Dirac equation. Recall that the leading spin structure 
of the f-jet is /j. We observe that the vertex 0^{l) is indeed generated by the 0{gs) expansion 
of $„(0, —oo)P±iIj{0), which is the part of the operator in Eq. ((7j) that sits to the left of the 
final state cut. 

To show that the operator in Eq. ((7j) gives the leading fragmentation contributions for Fl 
to all orders we proceed inductively by considering the coupling of two w-jet gluons to the hard 
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scattering, fig. 3. Because the hard scattering subgraph is UV dominated, its Dirac structure is 
again 7^ hence for Fl the same |^ factor enters at the current vertex. Then the graph of fig. 2 
is proportional to 



(T1T2) 0"H/l+/2) 



V ■ p V ■ p 

where terms that will vanish by the massless Dirac equation have been dropped. There is a 
hierarchy among the terms of the above expression once this factor is imbedded into a cross 
section diagram. The leading contribution comes from the first term and with its factor in 
square brackets eikonalized, i.e. from 



yX2 



{TiT2)0^'{h + k) y-9s^^j +(1^2) = 

-..^(T.T.) ((.,)-^ - + /.)-7]^^^) + (1 ^ 2). (15) 

All the other terms in Eq. (jl4j) are suppressed by at least one power of (1 — x) in the cross 
section, relative to the contribution that has been retained. We observe that Eq. (jl5|) is indeed 
the vertex in momentum space generated by the 0{g^) expansion of $t,(0, —oo)p±il){Q). More 
specifically, the first term corresponds to the piece that is 0{gs) in the Wilson line and 0{gs) 
in the covariant derivative and the second term corresponds to the piece that is 0{gg) in the 
Wilson line and 0{g^) in the covariant derivative. 

The above argument can be generalized straightforwardly for an arbitrary number of w-jet 
gluons attached to the hard scattering to the left of the final state cut. The leading contribution 
at a; — >• 1 is generated by ^^{0, —oo)p±'ip{0). This is how we show constructively, i.e. from 
analyzing the diagrams, that the new jet function J', as defined by the operator of Eq. ((7j), 
will contain the leading, 0{{1 — x)^), terms in the longitudinal cross section. Note that the 
operator in Eq. ((Tj) contributes at leading twist level. It is also the only operator that generates 
— contributions from the fragmentation of the final state jet. It can be shown that any 
other operator with higher number of covariant derivatives than in J' will genarate contributions 
of k = l,2,... 

Once the jet function J' is defined in terms of the expectation value of an operator it can 
be calculated in perturbation theory. However, the momentum integrals that are obtain in this 
calculation will also contain contributions in which the gluons coupling to the w-jet are soft. It 
is a well known (and non-trivial) statement that such contributions can be expressed in terms 
of expectation values of Wilson lines. Such contributions are removed by the denominator in 
Eq. ((7j) so that J' is by definition free from soft contributions. The latter ones will be contained 
only in the function V, Eq. (j2)). In fixed order perturbation theory this amounts to performing 
IR subtractions from the expressions generated by the numerator in order to obtain a J' jet 
function that contains only enhancements from final state fragmentation. 

There is also one more subtlety that arises when dealing with the function J', which is again 
characteristic of F^. The UV renormalization of the effective operator that defines J' captures 
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enhancements that are due to the near colHnear fragmentation of the outgoing jet. These contri- 
butions in the full theory (i.e. before factorization) are generated by genuine collinear infrared 
singularities. Although these singularities cancel upon adding real and virtual contributions in 
a full calculation, their finite terms remain and these are captured by the factorization theo- 
rem. So, collinear enhancements in the full theory are associated with ultraviolet singlularities 
in the effective jet operator and its ensuing anomalous dimension. In addition to these singu- 
larities, though, there are genuine ultraviolet singularities in the full theory that lead to the 
renormalization of its parameters, i.e. the coupling constant. Such ultraviolet singularies can 
also be present in the effective operator and these are not related to collinear enhancements 
but to the renormalization of the coupling constant in the jet function. If one is not careful to 
subtract these genuine ultraviolet singularities from the renormalization of the operator then 
one is in danger of double counting them when the renormalization group analysis is performed. 
In other words, the UV singularities of J' generate both the collinear enhancements and the 
renormalization of the effective coupling of the jet. This effect is not seen in low orders in F2 
because the jet function J starts at C(q;°), whereas in F2 it starts at 0{as). We shall return 
to this issue when we construct explicitly the anomalous dimension of J'. 



3.3 The RG equation for the new jet function and exponentiation 



The renormalization group equation for the jet and soft functions as well as the exponentiation 
of the logarithmic enhancements near the elastic limit are both consequences of the factorization 
formula. Therefore, we rewrite it here by making explicit the arguments that the participating 
functions depend on. Call p the total momentum of the "U-jet and (1 — z) the longitudinal mo- 
mentum fraction in the v direction that the jet carries into the final state. From the kinematics, 
X < z < 1, and = {l — z)Q'^/x is the small invariant mass of the jet. The jet has been defined 
with respect to the v direction hence the large scale on which it depends is (p ■ vy. Then the 
Fl factorization formula reads 



Hr 



{p ■ vY {p ■ v)" 



1^' 



L 



1 dz 



J' 



(16) 



Recall that p/ is the total momentum of the final state and that all functions are renormafized 
at scale //^. The only singularities present are the collinear ones with respect to the initial v 
direction. In moment space the above convolution becomes a simple product: 



FL{N,Q\e) = 



Hr 



1^' 



IJ? 



{p ■ V 



(17) 



where the Mellin transformation has been defined as usual 

dxx''-'FL{x,Q^,e), 



' jp ■ v f 



FL{N,Q',e) = r 
Jo 



dzz^'-^J' 



[l-z)Q' ip-v) 



(18) 
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and similarly for V and (f). Note that we have not included the overall factor into the 
definition of J'. This function will contain only the logarithmic in N terms. Also note that 
the UV dominated hard scattering function contains no dependence on N. Its presence 

is necessary, however, to render the right-hand side of Eq. (fTTjl independent of /x^. 

Following the approach of Collins, Soper and Sterman, (for review and references see ref. 0) 
we can write down the two equations that the jet function satisfies. The first is the covariance 
equation 



d 



d\n{p ■ v) 



In J' 



{p ■ 



as(/i^ 



K 



asif/) +G 



'{pj_u) 



-, as(/i^ 



(19) 



Since Fl is invariant upon variation of the light cone directions with respect to which the jet 
has been defined, the jet variations K and G are cancelled by the corresponding variations of 
the soft function V and the hard function \Hi\^ respectively. The dependence of V on the 
light cone directions has been left implicit in Eq. (fT7j) . The quark distribution function cf) is 
by definition independent of the final state direction v so it does not play a role in Eq. ()19|) . 
The function K is defined at scale /x^ and depends only on the common arguments of J' and 
V . The sum K + G is renormalization group invariant. The renormalization group equation 
for K reads 

d . oA 1 



-.K 



|7x(as(/x^)), 



with the well known cusp anomalous dimension in MS scheme [TT] 

5 



(20) 



CfCa 



67 
36 



TT 

12 



1^ 



CpNf 



(21) 



The second equation that the jet function satisfies is its own renormalization group equation 



d 



dlnjj' 



In J' 



{p ■ vf 



2' 



(22) 



We employ MS scheme which is a massless scheme, so the anomalous dimensions 'Jk and 7j' 
can depend only on the renormalized coupling. 
The solution of the covariance equation (|T9|) is 



In^ 



where we have used Eq. fj20|l and have defined 

= 7K(tts) + 2/3(a 



(23) 



B{as)/2 



d 



da. 

K{l,as)+G{l,as). 



with 



d 

din fj,' 



-a. 



47r 



A-kJ 



P2 + 



(24) 
(25) 

(26) 



10 



The solution of both the covariance and the RG equation for the jet J' is 



J'i 
In^ 



J'(l,l,a,((p-t;)2)) 



r.2 



r,2 



■7j'(«^(/^ )) 



and the boundary condition for the jet function is (setting {p ■ v)^ = Q"^ 

J'il, 1, as{Q')) = Cf^^ + 0{al{Q')). 

vr 

The renormahzation group invariance of Fl imposes the condition 

1H + 7J' + 7f + 70 = 0, 



(27) 



(28) 



(29) 



where the anomalous dimensions of all the functions involved are defined just as in Eq. (j22j) . The 
above relation allows us to write F^, Eq. ()17p . in a variety of exponentiated forms depending on 
which anomalous dimensions we wish to include in the exponent and what boundary conditions 
we choose for the factors Hl, J', V and 0. The form that is suitable for our purposes is the 
one that involves the jet anomalous dimension jji. Explicitly, 

FUN,Q',e) = ^\HL{l,l,a,{Q'))\\j\l,l,as{Q')){V-4>){l/N,a,{Q'),e) 



N 
X exp 

X exp 



dfi'^ ( Q 



d^l_ 

2/N /i^ 



In 



-A(«,(/i2)) + fi(a,(/i2)) 



d 



lj'{as{fi')) + 2[3{a,{fi'))—\nJ'{l,l,as{fi')) 

da. 



+ 0(-ln°iV 



(30) 



We shall be referring to this equation as the Sudakov exponentiated form of F^. By inspecting 
the above expression it is easy to see where the logarithmic in N corrections come from and 
how they are organized. As usual, the leading InA^ terms arise from the soft regions and are 
contained in the function A. The logarithmic corrections due to the fragmentation of the final 
state jet are contained in 7j/ and the ones due to the renormahzation of the coupling constant in 
the jet function are in beta function term. The beta function term arises from the evolution of 
the boundary condition from the low scale J'(l, 1, as{Q^/N)) to the high scale J'(l, 1, as(Q^)). 
It is worth empasizing again here that in Fl even the lowest order contribution cannot be 
normalized unambiguously because it starts at 0{as)- Since the renormalizability of the theory 
determines the structure of the beta function term in the exponent, we can absorb it into the 
anomalous dimension 7j/ as 



d 

7j'(«s) =7J'(as) + 2/5(a,)— lnJ'(l,l,a,) 

da.. 



(31) 



while keeping the boundary condition at the high scale, i.e. J'(l, 1, as(Q^)). The practical 
advantage of this simple redefinition is that when we compute the anomalous dimension of the 
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operator in Eq. (0) we do not have to subtract out the UV poles that renormahze the couphng 
as was suggested at the end of the previous section. In this fashion we compute 7j/ and not 7j/. 
Once the distinction between the two is made as in Eq. (jHTj) . there is no danger of overcounting. 
This is the approach that we follow in the next section. 

Not all logarithmic in dependence is contained in the exponential of Eq. pO|l . but it resides 
also in the V ■ factor. This can be captured too, if the corresponding anomalous dimensions 
are introduced in the exponent. We shall leave this factor as is, because it is common to Fl 
and F2 and contains no further information that is characteristic of the longitudinal structure 
function only. Our final objective is to connect the exponentiated form of Fl, Eq. ()30|) . with the 
operator product expansion factorization as given by Eq. (jH)). In the MS scheme the coefficient 
function in moment space and at scale /i^ = is given by 

Cl{N, QV/"' = 1) = ^-^^'^'''^ (32) 

f{N,Q\e) 

The regularized quark distribution function f{N,Q'^,e) can also be written as an exponential 
by simply integrating its DGLAP evolution equation in Q^, given one boundary condition. This 
is done in the presence of the regulator, i.e. in Z) = 4 — 2e dimensions. It is useful to choose 
for boundary condition the normalization [l] 

f{N,Q' = 0,e) = l, (33) 

which yields the solution 



/(iV,Q^e)=exp 



(34) 



Here Tgg is the Mellin transformation of the (singular) quark splitting function Pqq{x, e, as(/U^))- 
It is important to keep in mind that the expression in Eq. ()34|1 is defined only in the presence 
of the regulator e. The integration over the coupling constant down to /i^ = is also defined 
by using a dimensionally continued 

When we attempt to compute the coefficient function Cl from the ratio in Eq. (jH^ we 
encounter the following two problems. First, we cannot use the exponentiated form of Eq. (|3(J|) 
because Sudakov factorization organizes the exponent of F^ in terms of logarithmic in N de- 
pendence whereas the DGLAP evolution organizes the exponent of / in terms of logarithms 
of Q^. In other words the ratio Fl/ f is not a simple ratio of two exponentials because parts 
of the prefactors in the Sudakov exponentiation of Fl are contained in the DGLAP exponent 
of /. The consistent way to compute the coefficient function is to write also Fl as a single 
exponential by solving its own DGLAP equation in D = 4 — 2e dimensions. This is indeed what 
is done for F2. The boundary condition used in the F2 DGLAP equation is the normalization 

F2{N,Q' = 0,e) = l, (35) 

and the solution is again a single exponential, like in Eq. ()34j) . Then the coefficient function C2 
follows simply from Eq. ()32|) with L ^ 2. It is here that we encounter the second problem for 
the longitudinal strucutre function. Fl cannot be normalized as F2 above, because it starts at 
0(«.) 

^If we were to use the analytically continued ag that would give Fl{N, = 0, e) = 0. 
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Although we cannot treat formally Fl like F2, we observe that a large part in their respective 
Sudakov exponentiated forms is the same. Indeed, if we repeat the previous renormalization 
group analysis for F2, we get an answer for F2 which can be obtained from Eq. ()30|) by simply 
substituting F^, — > F2, H2 and J — > J' . F2 can be normalized as in Eq. (jH^ because 
J = 1 + 0{as)- Therefore, if we know the coefficient function C2 to some order in as we can 
predict the coefficient function Cl to one order higher by using the Sudakov factorization. All 
we need to do is to isolate the parts that are different between F2 and Fl in their exponentiated 
forms. The problem of normalization of F^ at = is solved by defining a function JF^ as 



Since the jet has been moded out, we can normalize the ratio as = 0, e) = 1. Then we 

compare the Sudakov exponentiation of F^ with the corresponding one for F2 and we identify 
the parts in the exponent that are different. These parts are the corrections that must be 
included in C2 to obtain Cl. Once we do this we restore F^ via Eq. (j36|) . The net result of this 
manipulation is that the coefficient function for Cl is predicted to be 

Cl = jr^2, (37) 

and the ratio of the two structure functions can be readily obtained from their Sudakov ex- 
ponentiated forms. The net effect of the ratio J^l/F2 is to remove from C2 the contributions 
coming from the jet J of F2 and replace them with the contributions from the J' jet of Fl- To 
0{a1) the above relation yields 

(AT, gV/i^ = 1) = {^K In^ N + l^'iP In at) + O (1 ln° n) . (38) 



4 Explicit calculations to 0{al) 

In this section we prove via explicit calculations the Fl factorization formula Eq. © to (9(a^), 
or equivalently that Eq. (jHHj) holds. This is done by comparing our results with the coefficient 
functions C^^\x, 1), n = 1, 2 in Eqs. (ini[I21)- The new object to be computed is the jet function 
J' defined in Eq. (|7j), since both the soft function V and the quark distribution are the same 
as in F2 and can be taken from the literature. We show that the singular ln(l —x) fragmentation 
terms in Cl (x, 1) do arise from the UV renormalization of the operator defining J' and we 
construct its anomalous dimension 7j, to 0{as). 

The analysis of subsection 3.2 resulted in a set of Feynman rules for computing the logarith- 
mic enhancements of the new jet function J'. Specifically, the coupling of gluons from the w-jet 
to the external lines is given by the effective non-local vertex 0'*'(/) defined in Eq. ()10|). In that 
expression, though, the transverse and light cone polarizations have been separated, whereas 
we find it more convenient for our calculations to compute traces with covariant objects. To 
this end we define the covariant vertex 0^(/)cov as 

o\i) = o\iu - U - l^^vA p. o\iu = f^' - 7-/ /) • (1) 
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Then 



tr 



V)- 



tr 



(2) 



so we can use the covariant vertex for computing traces. Note that in terms of the operator 
definitions, Ocov is generated by the operator in the numerator of Eq. (|7j) after performing a 
gauge transformation in the hght cone direction. All functions entering the factorization formula 
have been defined in a gauge independent way and from now on we employ the Feynman gauge. 
Finally, in the MS scheme it is convenient for presenting intermediate results to define the ratio 
of scales 



S 



-IE 



47rx/i 



(3) 



4.1 The 0{as) contribution to J' 

To first order in the coupling has no logarithmic singularities, the answer being simply 



Cpx. 



The Fl factorization formula to this order reads 

.(1) 



(4) 



(5) 



We check that J'^^^ obtained from Eq. ((7j) comes with the correct normalization to reproduce 
the above lowest order result. There is only one diagram to compute, that of fig. 4. From 
Eqs. (121) and (jH} we have 



1-e 



2 

rri-e 



^^''\l-2ej r(l-26) 
CFxS'[l-2e{\n{l-x) 



1)] 



(6) 



We note that the whole contribution to F^^ comes from J'^^^ and no IR subtractions are needed 
to this order. For convenience we have included in J'^^^ all overall normalizing factors associated 
with Fl- This short calculation shows that J' as defined by Eq. ((7j) is correctly normalized. 



4.2 The 0{al) contribut ion to the jet function J' 

The non-trivial infrared behavior of F^ starts at We compute the jet function to 

this order, J'^'^\ and show that after summation over final state cuts and IR subtractions as 
required by the denominator of Eq. ((7j) the only surviving singularities for the operator diagrams 
are simple UV poles. In section 3.3 we argued that the J' operator can be multiplicatively 
renormalized in moment space as 

j:(iV,Q2,a,(/i2),e) = Zj,{a,{fi'),e)J'{N,Q\as{fi')). (7) 

In a massless scheme such as MS, the anomalous dimension jj, can be constructed from the 
renormalization constant Zji as 

d 

ij,{as{ii^)) = -2as-^ResZj,{as{fi'^),e). (8) 
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Equivalently, to the order we compute, and because of the relations 

In AT 



dxx^-Hn(l- x) 
Jo 



N 



+ 0{ln°N/N), 



dxx 



N-l 



In^(l-x) 



-(lniV + 7s)2 + 0(ln°iV/iV), 



(9) 
(10) 



we can directly generate the (9(ln(l — x)) terms for J' in a;-space. Their coefficient is the 
sign-inverted residue of the UV singular part of J'^^-* . 

There are seven diagrams that determine J'*-^-* and they are depicted in fig. 5. Since cer- 
tain subtleties occur in the calculation arising from the presence of both linear (eikonal) and 
quadratic denominators in the loop integrations we proceed by analyzing each diagram in turn. 
In the full calculation to this order there are fourteen diagrams for F^. Whenever a diagram 
in our effective approach is free from IR divergences we provide direct comparison with the 
corresponding one of the full calculation taken from ref . . At the practical level, the effective 
approach is vastly more economical in capturing the logarithmic enhancements. 

Diagram (a): 

(a) ^ 2.(4^.)= d \ I iPS.. I £L (U) 

where a factor of 2 has been included from left<->right equal contributions. The numerator 
factor is 



(a) 



9l 



tr 



2 1 



p ■ p 



1 

— tr 



and the loop integrals yields through the use of Eq. Q in the appendix 



(12) 



2e 



(Pi 



(a) 



(27r)^ {l-pfP 
Then the diagram becomes 



-tr [0,^„, hi, 
9s 



-i 1 



i 1 



(An) 



-2.Aip-h)—-- + 0{e').{13) 



^2 ^2 ^ 



1 Q' 



e 2'KX Jo p ■ q p ■ q 



(47r)2 e 



4^C^1 + O(e0) 



(14) 



By including the overall factor from Eq. ((2|) we obtain the contribution to the unrenormalized 

J'(2) 



^1 j: 

TX 



1(2) 

u{a) 



ix 



47r g2 



TT 



^ lClx- + 0{e' 



This leads to the anomalous dimension 

lJ'(a){as 

as well as to the contribution for the renormalized jet function 



TT 



(15) 



(16) 



/(2) 
(a) 



-C|,ln(l -x) + reg. 



(17) 
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Diagram (b): 



and the numerator factor is 

1 



iCfCa 



2^ 



4£ 



dPSo 



(b) 



(27r)^ {p-lY{l-k2yP' 



N, 



9l 



tr 



' p-ki\ 



lQ{p-ki)l^ 1 



p ■ I 



(19) 



In the second step only terms that lead to UV singularities have been retained. There is a 
subtlety involved in the derivation of the above expression. A term proportional to l^lp is to be 
replaced as l^lp — *■ {l/2)l'^r]rjp. Note the factor 1/2 instead of the usual 1/4. The latter case is 
valid for the UV part of a loop integral if the integrand contains only quadratic denominators. 
In our case, however, it also contains a linear (eikonal) denominator and one must be careful 
not to overcount the polarizations that contribute to the UV divergent part (2 and not 4). 
Then, through the use of Eq. in the appendix we obtain 



2„ 2 
s 



iCpCA 



- fi'' / dPS2 I6{p ■ ki 



(47r) 



1 + ln 



P ■ k2 
p-q 



- + Oie' 
e 



aiQ' (CfCa\ 1 



+ 0(e°). 



Similarly to diagram (a), the corresponding contribution to the anomalous dimension is 

/ /X C^sCa 

and to the renormalized jet function 



t'(2) 



-Ci.C^ln(l-x) +reg. 



(20) 



(21) 



(22) 



Eq. reproduces the result for graph (10) in 0. 
Diagram (c): 



2i{4:7[as 



CfCa 



{2tt)d {i-pY(i-k^yp'' 



(23) 



and the numerator factor can be approximated in the UV as 

1 



:tr 



9t 



]{p ■ k 



1 



p ■ I 



(24) 



The same remark about counting the polarizations that contribute to the UV singularities as 
in diagram (6) applies here. Then, through the use of Eq. of the appendix we obtain 



(c) 



-8a^ \ C'f - 



Tlx 



Cl- 



2 

CfCa 



fi'' / dPS2 (p ■ k. 



\ + 2^1n^-^^~ 



p ■ k2 



p-q 



3\ 1 



+ 0(e° 



(25) 
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which contributes to the anomalous dimension as 



7j'(c)' 



C(2)-^ 



and to the renormahzed jet function as 



7'(2) _ 9 



CfCa 



m-l) Hl-x)+Teg 



(26) 



(27) 



Eq. pTjl reproduces the result for graph (13) in ff\. 

So far, the three graphs we have computed are in one-to-one correspondence with the 
corresponding graphs of the full theory. They are finite in the full theory and become UV 
divergent in the effective theory. Now we proceed to calculate the rest of the graphs. These 
are not in one-to-one correspondence with graphs of the full theory, which is the default case in 
factorization. No more comparison for intermediate results will be possible, only for the final 
answer. 



Diagrams (d) and 
tions yield 



e): For these two diagrams no subtleties are involved. Standard manipula- 



7j'(d)' 



TC 



j: 



'(2) 



id) 



+ reg. 



and 



7j'(e)l 



a. 



^ (--C^ + -Nj 
TT V 6 3 



t'(2) 



1 
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CFCA + -CFNf)\n{l-x)+Teg. 



(28) 
(29) 

(30) 
(31) 



The last two diagrams (/) and (g) are the most interesting to this order because they have 
both UV and IR singularities. Recall that the jet function J' is by definition free from IR 
singularities and its UV renormalization generates the fragmentation ln(l — x) terms in the 
coefficient function. We must therefore show that the left over IR singularities in the graphs 
(/) and (g), after summing over their two and three particle cuts, are precisely removed by the 
soft function V. The remaining UV poles will contribute to the J' anomalous dimension. 

Diagram (/): 

The two particle cut of this diagram (/2) is 



(/2 



2(47ra, 



iCfC^ 



2^ 



4e 



dPSo 



d^^l 



if) 



ajQ' [CfCa 
Attx 
'I 



S\l-x)- 



X 




(32) 



We have separated UV from IR singularities by retaining the sign in front of the regulator. 
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The three particle cut (/s) can similarly be obtained through the use of the three massless 
particle phase space given in Eq. (jH)) of the appendix in terms of convenient coordinates for our 
case. The result is 



(/s) = 2(-2)(47ra 



2 / iCFCA\ 1 ,,4e 

V 2 



/rfP^3(2-4)(l-e)(p-fci) 



+ 2k2) 



3 2 

+ 



(33) 



Adding the two and three particle cuts we obtain for the full cut diagram 

«2Q2 .CfCa 



if) = (/2) + (/s) 



Attx 
'1 

X 



SHl-x) 



-2S' l + -e 



+ 2 ln(l - x)-^ + 0(e°) 



(34) 



We note the cancellation of the double IR poles which arise from soft gluon emission. Simple 
IR and UV poles and their combination do survive. 

Diagram (g): 

Similar manipulations as in diagram (/) yield the following results for the two and three particle 
cuts. 



{92) 



2(-0(47ra,)2 



CfCa 



4e 



dPSo 



N, 



(9) 



2'KX 



2 CfCa 



S'(l-x) 



1 1 

7 + 



1 1 
+ 7 ■ 



-e) 



(^3 



2(47r«,)2 



CfCa 



/i^^ / dPS^ 16(1 - e) 



27rx 



CfCa 



S''H\-x\ 



-2fL 



(27r)^ il-k^fPip-l) 
+ 0(e°) 

(p-fci)(p- (fci + fcg)) 



(35) 



(36) 



Adding the two and three particle cuts we obtain for the full cut diagram 

CfCa^ 



{9) = {92) + {gz) 



2tix ^ 

X 



S'{l-x 
1 



+ ln(l - x) 



+ 0(e° 



(37) 



As in diagram (/), we note the cancellation of the IR double poles. UV and IR simple poles 
do remain. 

Finally, we add the contributions from diagrams (/) and {g) and multiply the result by 
{2'Kx'^ / alQ'^) ^ see Eqs. EI), to obtain their contribution to the unrenormalized J'^^\ It is 
convenient to separate the color structures, since they behave quite differently. The result is 



J, 



'(2) 



ClxS'{l-x) 



i - 1 + 



^ +ln(l-x)— ^ + 0(e°) 



+ 



CfCa 



xSHl-x) 



2e 
18 



+ 0(e° 



(3J 



We observe the absence of any IR poles in the CfCa color structure. 
Total contribution to J'^^-*: 

To O(a^) the jet function J'^"^^ is obtained from the sum of the diagrams {a)-{g) after subtracting 
the purely soft contributions to the same order. The IR subtraction term for starts at O(a^) 
and it is J'(^) O 1/^^). Explicitly, 



J'(2)(a;) = 42)(^)_(j/(l)^01))(^) 



(39) 



where the subscript u denotes that the jet function is unrenormalized and IR unsubtracted. 
The soft function 1^(1 — x) is universal and familiar. To 0{as) it is defined in fig. 6. In the 
Feynman gauge it is 



\/(i)(l -x) = CpS' 



^ 5{l-x)+ ^ 



1 



— e) V 1 — x 



+ 



ln(l — x) 
1 — X 



+ 0{e) 



(40) 



The argument of V is the longitudinal momentum fraction that enters the w-jet. The first 
term in the square brackets is of UV origin and comes from the virtual contribution. Note the 
absence of IR double poles. We compute the convolution J'^^^ V"*^^-* using the integral lemmata 
(jUmSl) in the appendix. The result is 



(J'« ®\/«)(x) 



J'^^\x)0V^^\l-x] 
ClS^' [x - ex(ln(l - 
1 



x] 



1)] 



-e e 



-6(1 -x) + 



ClS'^'x 



-(ln(l 

e 



x] 



X 



1) 



+ 



/ln(l 



XI 



\ 1 



X 



+ 0{e) 



i-e) 



ln(l -x) + C(e°) 



■ (41) 



As only diagrams (/) and (g) have non-trivial IR singularity structure, we consider only the 
(/) + id) ps^it of J'^^-* and perform the IR subtraction there. Specifically, 



j: 



/(2) 

{f+9)n 



X 



31 
'27 



CfCa 1 
2 



+ 0(e° 



(42) 



All IR singularities have been removed from J'*^^) as dictated by its definition in Eq. ((Tj). The 
left over terms are simple UV poles. Now we are finally able to compute the contribution of the 
diagrams (/) and {g) to the jet anomalous dimension and to the UV renormalized jet function 
in X space 



and 



t'(2) 



3 



( 3Cp 



TT 



2 



Cp + —CfCa 



ln(l — x) + reg. 



(43) 
(44) 



1 

.2^4 

It is true that we could have arrived at the above equation by just computing the UV simple 
poles of only the two particle cut graphs and converting them into ln(l — x) contributions. We, 
nevertheless, took the pain of calculating also the three particle cuts and the subtraction terms 
to show that indeed only simple UV poles survive in the final expression, thus showing that an 
anomalous dimension can be defined to sum the singular terms. 
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4.3 The final 0{al) result for Fl 

First, the total contribution to J'^^-* is obtained by summing the partial results in Eqs. ()17[ 
EHElEIlIll- It is 



J'(2) 



- - 2C(2) (C| - ^ 



6 



':CFNf 



ln(l — x) + reg. 



(45) 



The leading soft double logarithm is generated by the UV renormalized J'^^-* ® V^^\ Eq. ()4H1 . 
It is 

(J'(^) ® = ln2(l - x) +reg. (46) 

We have organized the summation of the singular ln(l — x) terms in such a way that all 
fragmentation (9(ln(l — x)) terms come from J'*^^^ and the soft (9(ln^(l — x)) term comes 
from J'(^) ® V^^\ At this stage we can determine the coefficient C^^\x, 1) as given by the Fl 
factorization formula. It is 



(x, 1) = (J'(^) ® \/(i))(x) + J"^^\x), 



(47) 



which, through Eqs. (j451 EUj) . exactly reproduces the singular part of the result from the full 
calculation in Eq. ()12|) . We have thus shown explicitly how the factorization formalism captures 
the logarithmic enhancements to (9(a^). Finally, the jet anomalous dimension 7j, is 



TT 



9 

—Cp — 2C A 
2 



4C(2)(C^-%)-^/3; 



(4J 



and the anomalous dimension 7j/ that does not include the running of the coupling is obtained 
from Eq. dM}. It is 



7j/(as) = — 

TT 



9 



-Cf - 2Ca - 4C(2) [Cf-^ 



(49) 



5 Summary and conclusions 

We have studied in detail the longitudinal structure function in deeply inelastic scattering near 
the elastic limit x — > 1 and have shown that it can be Sudakov factorized in terms of two jet 
and one soft functions. The soft function and the target jet function are the same as in F2, 
whereas the current jet function is characteristic of and measures the transverse momentum 
of the struck parton. We have shown that the quasi-elastic logarithmic corrections from the 
fragmentation of the current jet can be resummed in moment space in terms of the anomalous 
dimension 7j/, which we calculated to 0{as). 

We emphasize that our analysis is at the leading twist level but goes beyond the light-cone 
expansion. It resums terms that are of order In'^ N/N, k > 0. This approach can be extended to 
the case of F2, where such terms are present, but subleading relative to the order In*^ N, k >0 
terms. It is quite remarkable that terms that are power suppressed in N are amenable to a 
similar systematic analysis as the leading ones, just through the introduction of a single new 
function, the jet J'. This leads us to the conjecture that, indeed, all power in suppressed 
terms can be resummed in quasi-elastic processes, via the introduction of further jet functions 
defined by operators with higher powers of the transverse derivative. In Ref. a systematic 
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procedure has been presented for contructing these higher dimensional jet functions, and the 
present analysis is really a concrete example of that procedure. 

Returning to F^, we note also that the logarithmic enhancements from the final state 
fragmentation are delayed by one order of relative to the ones in F2. As these enhancements 
are summed to all orders, part of their effect (the beta function term in Eq. (j30|) ) is to change 
the overall normalization by shifting the scale of as in the preexponent factor from to Q'^/N. 
Thus, in the region x ^ 1 higher order corrections are parametrically as important as lower 
order ones. This has implications for obtaining the power corrections to F^, using for instance 
renormalon methods, since the power corrections come predominantly from this near elastic 
region. This enhancement, which is characteristic of F^, can be thought of as an indication 
from perturbation theory of the large power corrections that this observable is known to receive 
[TK] . These and applications of the formalism to other inclusive processes are currently under 
investigation. 
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Appendix 

• Three useful integrals with both quadratic and linear denominators, valid at = 

f d^l 1 -i 1 1 (p-k\ Q 

J (47r)^ (/ - p)2(/ - k)^p ■ I ~ (47r)2 e p ■ {p - k) [p ■ p ) ^ ^ ' 

0(e°) . (2) 



{2tt)^ {I ~ P p ■ I {A-nftp-p 



d^l _ 1 Vf.u . p-k 



ln^^ + 0(e°). (3) 



(47r)^(/-p)2(/-A;)2/2p./ {Att)^ 2e p ■ {p - k) p-p 
• The three massless particle phase space: 

The three particle phase space cannot be written conveniently in terms of the light cone mo- 
menta of the outgoing particles alone. In this case we prefer the parametrization in terms 
of Lorenz invariants defined in the rest frame of two of the three particles, say k2 and k^. 
Specifically, 

JO JO 

where 6 and (p are the polar and azimuthal angles of ki with respect to the p direction. In 
the ^2, ^3 CM frame, p defines the direction of the light cone, which is different from the one 
defined in Eq. (@]). The definitions of i/i and z are 

p-ki 0--x)-yi-{q-ki)/{p-q) 

yi = — zi = . (5) 

p-q {l-x){l-yi) 
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It follows that 

p-k2 , (l-cose\ p-ks /l + cose\ ki-k^ 

^ = (l-|/i)(^^^J, ^ = (i-yi)^^^J, — ^ = (6) 

For azimuthically symmetric integrands the d4> integration can be readily performed and in 
terms of the variable 

1 — cos^ 

-=—r- p) 

the phase space becomes 

X /' dyi y^il - yif-^' t dz, z^l - z,)-\ (8) 
Jo Jo 

Since all three particles are massless, the phase space is symmetric in their variables. So the 
phase space integrals in the ki, k^ CM system are obtained from the above expressions by 
permutation of the particle indices. 

• Finally, certain useful convolution lemmata: 

h{x)®5{l-x) = h{x). (9) 

= X (^ln(l -x) -lnx + i - 1^ . (10) 
x\n{l - x) ® (—^\ =a;ln2(l-x) + (l-x)ln(l- x) - -<(2). (11) 

\ -L X / -f- 

\ 1 — X J ^ 2 X 

X ® {1 + x^) {^^Y^^^ =xln^(l-x)+reg. (13) 
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Fig. 1. (a) Definition of the jet function J'. Dividing out the soft function V in rc-space is 
meant in the sense of inverse convolution (or convolutive division), (b) Fi factorization to 
all orders. 
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Fig. 2. The factorization of the w-collinear gluon (/) from the hard subgraph generates the 
effective vertex ®. Parallel lines around the photon denote contraction with p^. 





Fig. 3. The factorization of two t;-collinear gluons from the hard subgraph generates the 0{a 
expansion of ^yp±tp- Parallel lines around the photon denote contraction with p^. 
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Fig. 6. The 0{as) soft function V{1 — x) in the Feynman gauge. 
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